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THERMODYNAMICAL CONSISTENT MODELING AND ANALYSIS OF 
NEMATIC LIQUID CRYSTAL FLOWS 

MATTHIAS HIEBER AND JAN PRESS 


Abstract. The general Ericksen-Leslie model for the flow of nematic liquid crystals is reconsidered 
in the non-isothermal case aiming for thermodynamically consistent models. The non-isothermal 
simplifled model is then investigated analytically. A fairly complete dynamic theory is developed by 
analyzing these systems as quasilinear parabolic evolution equations in an Lp — L^-setting. First, 
the existence of a unique, local strong solution is proved. It is then shown that this solution extends 
to a global strong solution provided the initial data are close to an equilibrium or the solution is 
eventually bounded in the natural norm of the underlying state space. In these cases the solution 
converges exponentially to an equilibrium in the natural state manifold. 


1. Introduction 

The continuum theory of liquid crystals was developed by Ericksen and Leslie during the 1960’s in 
their pioneering work 0121]. This theory models nematic liquid crystal flow from a hydrodynamical 
point of view and reduces to the Oseen-Frank theory in the static case, see mm- It describes the 
evolution of the complete system under the influence of the velocity u of the fluid and the orientation 
configuration d of rod-like liquid crystals. Hence, d = d{t,x) is a unit vector in The original 
derivation Bm is based on the conservation laws for mass and linear as well as angular momentums. 
General liquid crystal materials are described by the Landau-de Gennes theory 0 from a unified point 
of view. 

The Ericksen-Leslie theory is nowadays widely used as a model for the flow of nematic liquid 
crystals, see for example the works of Ericksen and Kinderlehrer uni, Chandrasekhar, 0, deGennes 
and Prost 0 as well as Virga |40j . 

Note that these models are mostly formulated in an isothermal environment and are, in general, 
neither thermodynamically consistent nor thermodynamically stable. To the best of our our knowl¬ 
edge, only very few articles are dealing so far with the thermodynamical consistency of these models. 
Concerning the Ericksen-Leslie model, for a physically rigorous derivation we refer to the work of 
Muller [ 35 ] concentrating on the modeling aspect, and for recent analytical work to Feireisl, Rocca, 
Schiperna m, Feireisl, Fremont, Rocca, Schiperna m and Li, Xin |2^. Non-isothermal Landau-De 
Gennes nematic liquid crystal flows were investigated in the recent articles [14] and m- 

The aim of this paper is twofold: first, we reconsider the Ericksen-Leslie approach from the per¬ 
spective of thermodynamical consistency and stability. Following arguments from thermodynamics 
and employing entropy principles, we derive consistent models in a mathematically efhcient way, even 
in the case of compressible fluids. Let us emphasize that, in the end, our model contains the classical 
Ericksen-Leslie model in its general form as a special case. 

Secondly, we investigate our model analytically. Restricting ourselves to the case of constant 
density and not taking into account so called stretching, we develop a rather complete dynamic theory 
for the equations representing these models. More precisely, we first prove the existence of a unique, 
local strong solution to this system. We further show that this solution extends to a global, strong 
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solution, provided the initial data are close to an equilibrium or the solution is eventually bounded in 
the natural norm of the underlying state space. In this case the solution converges exponentially to 
an equilibrium in the natural state manifold. The results obtained thus parallel those proved recently 
by Hieber, Nesensohn, Priiss and Schade in [T^ dealing with the isothermal situation. For results 
concerning the asymptotic behaviour of solutions in the situation of the whole space we refer to 
the work of Dai and Schonbek |4] . 

The nowadays called simplified Ericksen-Leslie model in the isothermal situation was introduced 
and investigated first by Lin in |26| and m- Lin and Liu [55], [5S] studied the situation, where the 
nonlinearity in the equation for the director d is replaced by a Ginzburg-Landau energy functional. 
The existence of global weak solutions to this system in dimension 2 or 3 was proved under suitable 
assumptions on the intial data. For related results see [20] and m- Wang proved in [41] global 
well-posedness for the simplified system for initial data being small in BMO~^ x BMO in the case 
of a whole space by combining techniques of Koch and Tataru with methods from harmonic maps to 
certain Riemannian manifolds. 

The general Ericksen-Leslie model (in the isothermal situation) is based on the Oseen-Frank energy 
density functional which takes into account stretching as well as rotational effects for the director 
held. In the special case of homogeneous isotropic elasticity the equation for d reads as 

(1.1) dtd + u-Vd-Vd+^Dd=-^{M+\Vdld) + ^{Dd-d)d in(0,r)xll. 

Ai Ai Ai 

Here D = ^{S/u+[Vu]'^) denotes the symmetric, V = ^(Vu — [Vit]^) the anti-symmetric part of the 
deformation tensor and Ai,A 2 G R\{0} are material coefficients. Modihcations of this model were 
considered by Coutard and Shkoller in |3] in which the above equation (HU) for d is replaced by a 
Ginzburg-Landau type approximation. 

(1.2) ^{dtd + u ■ Vd + d ■ Vu) = Ad —^(|fi |2 — l)d in(0,T)xH. 

They proved local wellposedness for this system as well as a global existence result for small data 
within this setting. Note, however, that in this case the presence of the stretching term d ■ Vu causes 
loss of total energy balance and, moreover, the condition \d \2 = 1 in (0,T) x H, is not preserved 
anymore. For recent results on the general Ericksen-Leslie model with d satisfying (HU), we refer 
to the articles [Ml, m, Ha, m, isa, m and m], which contain well-posedness criteria for the 
general system under various assumptions on the Leslie coefficients. For results on non-isothermal 
models including the above mentioned stretching term, see m and references therein. 

Let us stress at this point that an important novelty of our approach lies in the fact that the 
complete model described in Section 2.7 is rigoroulsy proven to be thermodynamically consistent 
and stable. Specialising to the isothermal situation, we rediscover in particular the classical general 
Ericksen-Leslie system. It is interesting to compare our approach with the approach of Muller |33| . 
and with the energy variational approach developed by Liu and coworkers m and by Virga |40| . 

The plan for this contribution is as follows: Section 2 is devoted to the modeling of liquid crystals. 
In particular, based on the entropy principle, we derive a model of Ericksen-Leslie type which is 
thermodynamically consistent and stable. In Section 3 the equilibria of the system are identified - 
which are zero velocities and constant temperature and director - and it is proved that these are 
thermodynamically stable. The negative total entropy is shown to be strict Ljapunov functional, in 
particular the model is thermodynamically consistent. In Section 4 we prove local well-posedness of 
the non-isothermal simplified model and construct the resulting local semiflow in the natural state 
manifold of the system. We show that each solution which does not develop singularities in a sense 
to be specified converges to a unique equilibrium. These results are proved by means of techniques 
involving maximal Lp-regularity and results on quasilinear paprabolic evolution equations. For these 
methods, we refer to the booklet by Denk, Hieber, and Priiss |7] and to the work of Priiss and 
Simonett EH, Kdhne, Priiss and Wilke [l^, and LeGrone, Priiss and Wilke |23| . 
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By means of these techniques we are also able to prove analogous results for the full model, which, 
however, due to limitation of space will be presented elsewhere. 

2. Thermodynamical Consistent Modeling 

In this section we aim to give a self-contained presentation of a thermodynamically consistent modeling 
of liquid crystals. Like this we are able to refrain from refering to the orginal papers [9], mi, which 
are not easily accessible to a mathematical audience. Our approach does not only extend the classical 
Ericksen-Leslie model to the non-isothermal situation in a thermodynamical consistent and stable 
way but it also allows to exhibit the physical and mathematical beauty of this model. 

In this section, O C M" always denotes a domain with C^-boundary. 

1. First Principles 

We begin with the balance laws of mass, momentum, and energy. They read as 


dtp -I- div(pit) = 0 in O, 

(2.1) p{dt + u ■ V)u -\-VTT = d\v S in O, 

p(9t -I- It • V)e -I- divg = S' : Vu — Trdivu in O, 

u = 0, q ■ ly = 0 on dil. 


Here p means density, u velocity, tt pressure, e internal energy, S extra stress and q heat flux. This 
immediately gives conservation of the total energy. In fact, we have 

p{dt + u ■ V)e -I- div{q + ttu — Su) =0 in fl, 

where e := |mP/2 -I- e means the total mass specific energy density (kinetic and internal). The energy 
flux is given by := g -I- ttu — Su. Integrating over H yields 

i9tE(t) = 0, E{t) = Ekin{t) + Eint{t) = [ p{t,x)e{t,x)dx, 

JQ 

provided 

(2.2) q ■ u = u = 0 on dEl. 

Hence, if (12.21) holds, total energy is preserved, independent of the particular choice of S and q. 

2. Thermodynamics 

Assume a given free energy if of the form if = 'if{p, 9, r), where 9 denotes the (absolute) temperature 
and r will be specified later. We then have the following thermodynamical relations: 

€ = + 9r] internal energy, 

(2.3) p = -defy entropy, 

K = dge = —9d‘g'if heat capacity 

Later on, for well-posedness of the heat problem, we require «: > 0, i.e. '0 to be strictly concave with 
respect to 0 S (0,oo). 

In the classical case, where if depends only on p and 0, we have the Clausius-Duhem equation 
p{dt -E u ■ S/)p -E div(g/6() = S : yu/9 — q ■ V9/9^ -E {p^dpif — 7r)(divu)/0 in Q. 

Hence, in this case the entropy flux <I>^ is given by <I>^ := q/9 and the entropy production by 

9r := S : Vu — q ■ V9/9 -E {p^dpif — 7r)(divu). 

Employing the boundary conditions (12.211 , we obtain for the total entropy N by integration over El 

dtN{t) = / r{t,x)dx>0, N(t) = / p{t,x)p{t,x)dx, 

Jn Jn 


4 


MATTHIAS HIEBER AND JAN TRUSS 


provided r > 0 in SI. As divrt has no sign we require 

(2.4) TT = p'^dp'ip, 

which is the famous Maxwell relation. Further, as S and q are independent, this requirement leads 
to the classical conditions 

(2.5) S' : Vu > 0 and 9 • V 6 » < 0. 

Summarizing, we see that whatever one chooses for S and q, one always has conservation of energy and 
the total entropy is non-decreasing provided (12.511 . (j2.411 and (j2.2|l are satisfied. Thus, these conditions 
ensure the thermodynamic consistency of the model. 

As an example for S and q consider the classical laws due to Newton and Fourier which are given 

by 

S := Sn ■= “ipLsD + ^hdivu/, 2D = (Vu -I- [Vu]^), q = — aoV0. 

In this case, (12.511 is satisfied as soon as pig > 0, 2/is -|- npb > 0 and ao > 0 hold. Note that it does not 
matter at all whether fj,s,p,b,Oio s^re constants or whether they depend on p, 0, or on other variables. 

3. Nematic Liquid Crystals 

For isotropic nematic liquid crystals we assume a free energy density '0 of the form 

V'=' 0 (p, 6 »,t), with r = i|Vd|^. 

Here d means the orientation vector, also called the director, which should satisfy the condition 

n 

\d\l = 

Note that 2t = tr(V(i[Vd]^) is the first invariant of the matrix Vd[V(i]^, and for its last invariant it 
holds det(Vd[Vd]^) = (detVd)^ = 0, as Vdd = 0 by \d \2 = 1. 

We neglect spin energy below but take into account transport of energy due to couple stress. This 
means that the energy flux is replaced by 

$e '.= q + TTU — Su — IVDtd, Vt = dt + u ■ Vd, 

where H has to be modeled. 

As constitutive laws we will employ 

(2.6) S = Sn + Se + Sl, =-AVd[Vd]'^, q = -a^Ve - ai{d ■Ve)d. 

Sn means the Newton stress introduced above, Se the Ericksen stress, and Sl the Leslie stress which 
will be defined later. Assuming these two constitutive laws we derive in the following the balance of 
entropy, i.e. the Clausius-Duhem equation. A short computation gives 

(2.7) p{dt-\-u-V)!]div^rj = r, 
with = qjo, and 

Or = —q ■ VO/O + 2pLs\D\l + pb\divuf {p^dpij — 7 r)divu 
-k {pdrip - A)Vd[Vd]'^ : Vu -k (H - pdripVd) : Wtd 
+ Sl '■ Vu (divH -(- /3d) ■ Vtd. 

for some scalar function /3. Note that d • TAjd = 0 as |d |2 = 1, hence /3 S K. can be chosen arbitrarily. 
For the entropy production r to be nonnegative, we require 


Ps > 0, 2ps -I- npb >0, ao > 0, ao -I- ai > 0. 
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Except for the last one, these conditions are the well-known conditions from fluid dynamics, see 
Section 2.2. The subsequent terms in the definition of r have no sign, hence we require them to 
vanish, which yields the relations 

(2.8) TT = p^dpip, A = pdr^! n = pdripyd. 

Finally, to obtain nonnegativity of the last two terms, in the simplest case, we may assume that the 
Leslie stress Sl vanishes, and 

'^Vtd = div[(p9T-'0)V]d -f fdd, 

for some 7 = j{p,9,t) > 0. The condition |d |2 = 1 then requires /3 = AjVdp, which leads to the 
equation 

(2.9) j(dt +u-V)d = div[AV]d-f X\Vd\^d, 

a nonlinear convection-diffusion equation for d. This is the basic equation governing the evolution of 
the director field d. With these assumptions the entropy production reads as 

Or = -q- VO/9 + 2ps\D\\ + Pbldivup -I- 

7 

where 

a = div[AV]d -)- X\Vd\‘^d = 'jVtd. 

At the boundary did, energy should be preserved, which means ihe • u = 0. As g • u = 0 and u = 0 
this yields 

A9j/(i • dtd = 0. 

This is clearly valid if d satisfies the Neumann condition d^d = 0, which is physically reasonable. 

4. Stretching and Vorticity 

Observe that the equation (12.91) for d admits the solutions d = const, no matter how the velocity field 
and the temperature field are defined. In this case the director field is not at all affected by the fluid 
dynamics. This seems to be physically unrealistic and so the model should be adapted. 

This can be done by introducing a so-called stretching stress. To introduce this stress we follow 
Leslie. Define = I — d ® d the orthogonal projection onto := {d}'*', the vorticity V according 
to 2V = Vu — [Vw]^, and set 

n = pvVd + poPdDd — ^Vtd, 

where pv, are scalar functions of p, 0, r and 7 > 0. For brevity we use the notation 

a := Pddiv(AV)d = div(AV)d -I- A|Vdpd. 

Now we define the stretch tensor 

(2.10) gstretch^ 

27 27 

This modification of the model does not change the entropy flux = q/0, and the relevant entropy 
production becomes, 

gstretch :^u + Vtd-a = -(|n|2 jVtd ■ n)) + a ■ Vtd 

7 

= -(|ap -f- (n -7 a) ■ {pyVd + poPdDd - a). 

7 

If we want to keep the total entropy production at the same level as in the previous section, the 
simplest way to achieve this is to set n -7 a = 0, which yields the equation 

(2.11) "i{dtd -7 u ■ Vd) = div(AV)d -7 A|Vdpd -7 pyFd -7 poPdDd. 
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This is the stretched equation for d. Note that it preserves the constraint |(i |2 = 1- The entropy 
production is the same as before, we have 

9r = [ao|V 0|2 + ai{d\V9)‘^]/9 + 2fis\D\l + /Jbldivwp + -\Pddiv{W)d\l. 

7 

In particular, N satisfies 9tN(t) = f^r(t,x)dx, and so — N will shown below to be a strict Lyapunov 
functional for the system, as soon as 

(2.12) /is > 0, 2/is + nfib >0, ao > 0, ao + «! > 0, 7 > 0, 

and 

(2.13) K > 0, A > 0, dpTT> 0. 

Note that no conditions on the new parameter functions are needed, so far. 

5. Additional Dissipation 

We may add additional dissipative terms in the stress tensor of the form 

(2.14) = ^(n( 8 )d + d( 8 )n) + (PaDd ^d + d(^ PaDd) + MDd\d)d <» d, 

7 27 

where as before 2D = Vu + [Vu]^ and 217 = Vu — [Vm]^ are the symmetric and antisymmetric 
parts of the rate of strain tensor Vu. Note that the tensor is symmetric. Adding these terms 

to the stress tensor will be thermodynamically consistent provided their contribution to the entropy 
production ensures that the total entropy production remains nonnegative. By a simple calculation 
we obtain 

2 

(2.15) : Vu = : D = 2^{n\PdDd) + (/xl + f^)\PdDd\l + fioiDd\df. 

7 7 

So with n = — a, the total relevant dissipation amounts to 

-(|a |2 + 2fj,p{n\PdDd) + n^lPdDdH) + ^j,L\PdDd\l + no{Dd\df 
7 

= -|a - fipPdDd\l + fiL\PdDd\l + fio{Dd\df, 

7 

hence the total entropy production becomes 

9r = [ao|V 0|2 + ai{d\S/9)'^] / 9 + 2^s\D\2 + /ibidivup 

+ -|Pddiv(AV)d - ^ipPdDd\l + HL\PdDd\l + fio{Dd\d)'^. 

7 

Note that so far the parameter functions /x/, j = 0, s, b, V, D, P, L, ao, cJi, and 7 for thermodynamical 
consistency are only subject to the requirements 

(2.16) ao,ao + ai>0, /x^, 2/Xs + n/X{, > 0, /xo,/xl>0, 7 > 0. 

Recall that all parameters functions are allowed to be functions of p, 9, t. 

Remark, (i) A more refined algebra shows that it is enough to require 

2/Xs + pl > 0, 2/Xs + /xo 7 0 

in the incompressible case, and additionally 

do ^ to , , , Po. 

—2 ^ (2/Xs + Po)[ - \- Pb - n) 


in the compressible case. 
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(ii) We want to stress that in case fiv = 7 , our parameters /is,/xo,/lu,M d,/ ip,M l are in one-to-one 
correspondence to the famous Leslie parameters ai,..., uq. This shows that our model contains the 
isotropic Ericksen-Leslie model as a special case. 

6. Conservation of Angular Momentum in 3D 

We briefly discuss conservation of momentum in the physical important three-dimensional case. Recall 
that the mass specific density m of angular momentum is defined by 

m = a; X u. 

Balance of angular momentum reads as follows 

Vt{pm) + div(pu X m — a; X T) = — Cj x Tj, 

where we use Einstein’s sum convention and Ti denotes the i-th row of the stress tensor T. Thus the 
flux of angular momentum is given by 

= pii X m — a; X T. 

It is well-known that x = 0 in case T is symmetric. Therefore we may concentrate on the 
non-symmetric part of T which is given by 

^(nOd-dOn). 

27 

This implies 

Ci X = ^(n X d — d X n) = — — d x n 
27 7 

= — d X div(AV)d 

7 

= d^i^Xd X d,d) - 5,(^)Ad X a,d - ^d,d X d,d 
7 7 7 

= d^i^Xd X d,d) - d,i^)Xd X d^d. 

7 7 

This shows that e* x is a divergence provided = 0, i.e. if 

P-v = co 7 , for some constant cq € K. 

Then the flux of angular momentum becomes 

=pMXm—a;xr-|- cgXd x Vd. 

We mention that if we also require so-called objectivity of the model, then cq = 1, which means 
Pv = 7- 

7. The Complete Model: non-isothermal, compressible fluid, isotropic elasticity 


Summarizing, the complete model may be represented as 

dtp + div(pM) =0 in n, 

p{dt + u ■ V)u + Vtt = div S in 12, 

(2.17) p(5t -I- M • V)e -I- divg = S': Vu — TrdivM -I- div(AVdI>td) in 12, 

^{dt + u ■ V)d — pvVd — div[AV]d = A|Vdpd -I- pnPdDd, in fl, 

u = 0, g • u = 0, di,d = 0 on 512 . 

p{0) = po, u{0) = uo, 0(0) = 00 , d(0) = do in 12 
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These equations have to be supplemented by the thermodynamical laws 

e = %p + 0r], ?7 = -dgip, n = dee, 

(2.18) TT = p^dpTp, \ = pdri^/9, 
and by the constitutive laws 

S = Sn + Se + 

(2.19) Sn = 2psD + pbdiYuI, Se = -0XVd[yd]^, 

gstretch ^ d:D+J^^ PD - PV ^ ^ ^ ^ ^ poPdDd - jVtd, 

27 27 

gdiss ^^(n(g)d + d( 8 )n) + + I^P (^d + d^ PdDd) + po{Dd\d)d 0 d, 

7 27 

q = -aoVe -ai{d\Ve)d. 

Here all coefficients dj,oij and 7 are functions of p,0,T. For thermodynamic consistency we require 

( 2 . 20 ) Ps>0, 2ps + npb>0, ao > 0 , ao + ai>0, /xo,Ml> 0 , 7 > 0 , 

Finally, we will use in addition the following conditions 

Ps > 0 , 2ps + npb > 0 , ao > 0 , ao + ai > 0, 7 > 0 , 

(2.21) K > 0, A > 0, dpTT> 0, 

to identify the equilibria and to investigate their thermodynamic stability in Section 3. 

8 . The Complete Model: non-isothermal, compressible fluid, non-isotropic elasticity 

For the sake of completeness, we comment briefly on the non-isotropic case. Then ip = ip{p, 9, d, Vd), 
and the Ericksen stress tensor becomes Se = ■ Following the derivation in Sections 2.3 

and 2.4, here the energy and entropy fluxes read again as 

:= g -1- TTM — Su — UT>td, ^ri = qjd, 

and the equation for d becomes 

j'Dtd=Pda, a = diipS/did^ip) - pX/d'tp, 

in the case without stretching, and 

-fVtd = Pda -I- pvVd + poPdDd 

in the stretched case. The couple stress here is H = pd^dtp, and the entropy production now reads as 


9r = [oolVdll + ai(d|V0)^]/0 -I- 2ps\D\\ + /Xb|divup 

-h -|fyi(a - ppDd)\l + pL\PdDd\^ + po{Dd\d)‘^. 

7 

Summarizing, the complete model in the case of non-isotropic elasticity becomes 

dtp + div(pM) =0 in O, 

p{dt + u ■ V)u -I- Vtt = div S in O, 

(2.22) p{dt + u ■'\7)e + div q = S : Vu — Trdivu + div{pdydi’Ptd) in O, 

'y{dt + u ■ V)d - Pd(div(p^^) - pS/dip) = pvVd + poPdDd, in H, 

u = 0, q ■ p = 0, = 0 on dH. 

p(0) = po, u{0) = uo, 9{0) = 9o, d{0) = do in H 
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These equations have to be supplemented by the thermodynamical laws 

(2.23) e = ip + 9r], rj = -dgtp, k = dee, tt = p^dpip, 
and by the constitutive laws 

(2.24) Sn = ‘2nsD + ^jLeA\YuI, Se =-p-^-^\S/dY, 

gstretch ^ PD+PV ^ ^ ^ ^ PD - PV ^ ^ ^ ^ y ^ioPdDd - ^Vtd, 

27 27 

_ 1 _ 2 

gdiss ^^(n(8)d + d(8>n) + {PdDd ®d + d® PdDd) + no{Dd\d)d 0 d, 

1 27 

q = -aoVe - ai{d\Ve)d. 

Here all coefficients fdj,aj and 7 are functions of p,0,T. For thermodynamic consistency we require 
as before only ( 12 . 201 ) . We also note that the natural boundary condition at dVl here becomes 

ViVo^di’ = 0 . 

Observe that this condition is fully nonlinear, in general, in contrast to the isotropic case. 

Concluding, we mention as an example the classical Oseen-Frank free energy density for the isother¬ 
mal incompressible case, which is given by 

= ki{divd)^ + k 2 \d X {V X d )\2 -I- fcajd • (V x d)\^ + (^2 -|- fc 4 )[tr(Vd)^ — (divd)^], 
where ki are given constants. 

3. Thermodynamical Consistency and Stability 

In this section we determine the equilibria set of the complete system described above in Section 
2.7, show that the critical points of the entropy functional coincide with these equilibria and prove 
that they are thermodynamically stable. We begin investigating the set of equilibria. 

1. Equilibria 

Suppose that in some time interval t G (tijO) we have cltN(t) = 0. Then r > 0 implies r{t,x) = 0 in 
H. This yields V0(t, x) = 0 in H as ao > 0 and ao + ai> 0. Hence, 0 = 0* is constant in H. 

Next, by ps > 0, 2ps + npe > 0, we also have D = 0 in il. By Korn’s inequality and the no-slip 
boundary condition for u, we hence obtain u = u* = 0 in H, t € (tijO)- Therefore dtp = dtu = 0, 
which implies Vtt = 0 . 

Finally, 7 > 0 yields Htd = 0 in H, which implies that d satisfies the nonlinear eigenvalue problem 


div(a(a;)V)d -I- a(x)\'Vdl^d 

= 0 

in n, 

\d\2 

= 1 

in n, 

d^d 

= 0 

on dn, 


where a{x) = X{p{t,x),9{t),T{t,x)), for each fixed t G (^ 1 ,^ 2 )- But, as the next lemma shows, this 
implies Vd = 0 in H. Hence, d = d* is constant. 

Lemma 1. Let q > n, a G a > 0 and suppose that d G 77^(12;®.") satisfies (13.11) . Then d is 

constant in H. 

Proof. The idea is to reduce inductively the dimension N = n oi the vector d. This can be achieved 
by introducing polar coordinates according to 

di = Cl cosif, d 2 = Cl sin(/?, dj = cj-i, j > 3. 

Simple computations yield 

l = \d\l = \c\l \Vd\l = \Vc\l + cl\V^\l 
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and 


div(aV)cj + a[|Vc |2 + =0 in ^2, 

as well as dyCj = 0 on 912 for j = 2,..., n — 1. Moreover, by some more calculations we further obtain 
—div(aV)ci + acl\V^p\\ = a[|Vc |2 + Ci|V(/j|f]ci in f2, 


and 

as well as 


cidiv(aV)</5 + 2aVci • Wtp = 0 in f2, 


d^ci = cid^ip = 0 on 9r2. 

Multiplying the second of the last equations by ci(f and integrating over f2 we deduce 


0 = / [cidiv(aV)</3 + 2aVci • V(/3)]ci</3dT = / div[ciaV(p](pdx = — / Cia\V(p\^dx. 

J vl o j ^ 

Hence, CiV</5 = 0 as a > 0 by assumption. This implies that c satisfies equation m, where the 
vector c has dimension TV — 1. Inductively, we arrive at dimension = 1 and if d is a solution of 
mi with dimension 1, then d=lord=—Iby the connectedness of 12. □ 


Knowing that 9 and d are constant in 12, and Vtt = 0, we see that tt = p^dpipip, d, 0) is constant, 
hence p = p* is constant, provided the function p i-A 7 r(p, 0,0) is strictly increasing. This shows that 
we are at an equilibrium (p*, u*, 0*, d*) € S with 

^ = {(p*,M*,d*,d*) e (0,oo) X {0} X (0,cx)) X R” : |d *|2 = 1}, 

the set of physical equilibria. In particular, the functional — N is a strict Lyapunov functional. 

Observe that £ forms an n + 1-dimensional manifold. If we take into account conservation of mass 
and energy. 

Mo := / pdx = p*|I2|, Eq := / (p|u|2/2 -I- pe)dT = p*e*|I2|, 

Jn Jn 

at an equilibrium, then the values of p* and 0 * are uniquely determined by 

P* = Mo/|I2|, e* := e(p*,0*,O) = Eo/Mo, 
whenever 0 i-A e(p, 0 , 0 ) is strictly increasing, i.e. whenever k > 0 . 


2. Critical Points of Total Entropy 

a) Consider the entropy functional N with constraints of prescribed mass M = Mq and energy E = Eg, 
as well as G(d) := (|d|| —1)/2 = 0. Suppose we have a sufficiently smooth critical point (p, u, 0, d) of N 
with p,d > 0, subject to the constraints. Then the method of Lagrange multipliers yields km, ke G R. 
and KG S L 2 (£ 1 ) such that 

(N' + kmM' + keE' + kgG'\z) = 0, 

where 2 : = (ct, v, d, d). We have 

(M'|z) = / adx, {kgG'\z) = f Kcd-Sdx, 

J fl Jfl 

and 

(N'l^) = f [{dp{pr]))a + pdeipd + pdrpVd : VS]dx, 

Jn 

(E'|z) = f [pu ■ v{dp{pe))(T + pdged + pdreVd : VS]dx. 

Jn 


as well as 
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This yields the relation 

0 = / {[dpipni) + Km + ke{\\u\ 1 + dp[pe))\a +[pd0'q + KEpd0€]'d}dx 

Jn ^ 

+ / {kePU • V + [pdrp + KEpdr^lVd : S/S + Kad ■ 6}dx. 

Jn 

We first vary d to obtain p{dop + KEdge) = 0, which by p > 0 and by the definition of rj, e and k > 0 
yields ke = —IJO. Hence, 9 is constant and ke < 0. Next, varying v we obtain u = 0, as ke and p 
are not zero. Next we vary <5, which after an integration by parts, employing the boundary condition 
dud = 0, implies 

div(AV)d + Kcd = 0 in H. 

But then \d \2 = 1 implies kq = A|Vd||, and d is a solution of the problem (I3.1L which by Lemma 
1 shows that d is constant. Finally, we vary a to the result that dp{pip) = 9km is constant. As 
TT = p^dpij) is strictly increasing in the variable p, this shows that p is constant in H as well. Therefore, 
the critical points of the entropy functional are precisely the equilibria of the problem. 

b) Let 

H := N" + kbE" 

denote the second variation of N. Note that M" = 0 and kq = A|Vd |2 = 0. The identities 
Pidr-q - = -A, dpdg{pq) - ^dpde{pe) = 0, 

dl{pq) - \dl{pe) = d^q - ^d^e = 

imply 

-(ddzlz) = + A|Vd|^]da; > 0, 

by K, A, dpTT > 0. This shows that the second variation of N at an equilibrium is negative semi-definite, 
which means that the equilibria are thermodynamically stable. 

c) Summarizing we have the following basic result 

Theorem 1. The complete model has the following properties. 

i) Along smooth solutions total mass M and energy E are preserved. 

ii) Along smooth solutions the total entropy N is non-decreasing. 

Hi) The negative total entropy is a strict Lyapunov functional. 

iv) The condition \d \2 = 1 is preserved along smooth solutions. 

v) The equilibria are given by the set of constants 

£ = {(p*, 0, 9.f,dt,) : pt,,9.f e (0, oo), d* e R", |d *|2 = 1}. 

Here p*, d* are uniquely determined by the identities 

p* = Mo/|H|, e(p*,6>*,0) = Eq/Mq. 

vi) The equilibria are precisely the critical points of the total entropy with prescribed mass and energy. 

vii) The second variation of N with given mass and energy at equilibrium is negative semidefinite. 

In particular, the model is thermodynamically consistent and it is also thermodynamically 
stable. 

3. The Isothermal Case 

In the isothermal case we set 9 = const and ignore the equation for the energy. In this case, instead 
of the total mass specific energy e = 11112/2 + e, we employ the available energy Ca which is defined 
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by Ca = \uW/2 + ij}. We have the following balance of Ca which is a direct consequence of balance of 
total energy and entropy 

p{dt + u ■ V)ea + div($e — 0$,,) = —Or — prjVtO — • V0. 

In the case where 0 is constant this reduces to 

p{dt + u ■ V)ea + div($e - = -ra, 


with 


Ta = 2ps\D\l + ^t,|divu|^ + po{Dd\d)'^ + pL\PdDd\l + |Pd(div(AV)d - ppDd)\l/j. 

Therefore, in the isothermal case, the total available energy Ea is a strict Ljapunov functional for the 
system, i.e. 

dtEait) = - ra{t,x)dx, Ea(t) = / p{t,x)eait,x)dx. 

Jn Jn 

As a consequence, the equilibrium set is the same as in the non-isothermal case, dropping temperature, 
hence is a manifold of dimension n, and when we incorporate preserved mass it is isomorphic to the 
unit sphere in R". In this case the equations read 

dtp + div(/ 9 u) =0 in fl, 

(3.2) p(9t + u • V)m + Vtt = div5' in 11, 

"f{dt + u ■ y)d — pvVd — div(AV)d = AlVdIad + poPdDd in 


where ip = ip{p,T), tt = p^dpip, A = pdrip/0, and 


S = 2psD + pbdivul- 9XVd[Vd] ' + S-f+ S'; 


(3.3) 


(3.4) 


S' 


stretch dD + PV _ 


Pd — Pv 


27 


27 


d( 8 )n, n = pvVd + ppPdDd — "fVtd. 


S 


diss _ PP 


= 

7 


I d + d 0 n) + 


7PL + Pp 
27 


{PdDd ( 8 > d + d (g) PdDd) + pti{Dd\d)d ® d. 


If one further restricts to the incompressible case p = const > 0, A constant, /jq = po = Pv = Pp = 
Pl — 0, with p = ps one obtains the so-called isothermal simplified Ericksen-Leslie model 

p{dt + u ■ V)m -I- Vtt = /igAit — Adiv(Vd[Vd]^) in Q, 

(3.5) |d |2 = 1, divrt = 0 in fl, 

7(i9t-I-u • V)d — AAd = A|Vd|^d in fl. 

Of course, in all cases we have to add initial conditions as well as boundary conditions u = d^d = 0 
on dLl. Problem (j3.5l) subject to the condition |d |2 = 1 in O has been analyzed in a fairly complete 
manner in the recent article [18] by Hieber, Nesensohn, Priiss and Schade. 


4. Analysis of the Non-Isothermal Simplified Model 

In this section we consider the incompressible case p = const and we let p = ps- Hence, the 
pressure tt is no longer determined by Maxwell’s relation; it is now a free variable, a Lagrangian 
multiplier to cover the constraint div u = 0. Furthermore, in the following we neglect stretching, i.e. 
we assume pn = Pv = pp = Pl = Po = 0. For simplicity we also set oi = 0 and a = uq. Then the 
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resulting model - which we call the non-isothermal simplified Ericksen-Leslie model - reads as follows. 
pVtU — 2div(^Z3) + Vtt = div(AVd[Vd]^) in fl, 

|d |2 = 1 , div u = 0 in Sd, 


pKVfO - div(aV0) = 2p,\D\l + 6»AVd[Vd]'^ : D 
(4.1) - pdreVd : VtVd + div(AVdX't) 

'yVtd — div(AV)d = A|Vd||d 
u = d^O = dvd = 0 
■u(O) = Mo, 6»(0) = 6>o, d{0) = do 

Recall that p > 0 is constant and a, 7 , p as well as A = pdrf’/O^ k = dge 
0 > 0 and T > 0 . 


in LI, 
in LI, 
on dLl, 
in LI. 

—Odgip are functions of 


1. Regularity Assumptions (R) 

The parameter functions should have the following minimal regularity properties: 

/i, a, 7 e ^^((0, 00 ) X [0, 00 )), i/: € 00 ) X [0, 00 )); 

We also require the positivity conditions 

p > 0, a > 0, K > 0, 7 > 0, A > 0, 

which have been mentioned before, but for well-posedness of the problem for d we need in addition 
to require A + 2 t3t-A > 0. We assume that LI C M" is a bounded domain with C^^-boundary. 

2. Maximal Lp-Regularity of the Principal Linearization 

The equation for u will turn out to be only weakly coupled, so we first concentrate on the system for 
w := \9, d]^. The principal part of the linearization becomes 

dtwA{wo,'V)w = f inf], 

(4.2) di,w = 0 on dLl, 

w(0) = Wo in LI. 

The matrix A = A{wo, V) reads as 

, — bo^do : (AqA T (9.7-Ao[VdQ]'^VdQ : V^)V 

6i[Vdorv, -7o-i(AoA + 9,Ao[Vdo]^®Vdo:V2). ■ 

Here kq = k(0o,to) etc., and we used the abbreviations 

ao pdoldrPof , ddrPo , pdrPo 

ao = -, ai =-, bo = -, 0 i =-. 

PKo 7oKo 7oKo 7o 

Note that A(wo, V) is second order in the diagonal, but third and first order off-diagonal! This is a 
mixed-order problem subject to Neumann boundary conditions and subject to variable, non-smooth 
coefficients. For resolvent estimates within the L^-setting for various mixed-order systems we refer 
to the work of Grubb HZ]. Regarding the maximal -regularity of this nonstandard problem, we do 
not know of any general theory covering the above situation. However, we note for the whole space 
case there is the theory of Denk and Kaip [5] available. Nevertheless, we prove maximal Lp-regularity 
for this problem in the following. 

To this end, fix q S (1, cxd) and choose as a base space 

To := Lg{Ll)xHl{Ll-,W^), 

and as a regularity space 

Fi := {w = {9,d) e Hl{Lt) x H^Ll-W) : 8,^9 = = 0 on dLl}, 
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equipped with their natural norms. We will also employ the time-weighted spaces defined by 

m G No, p G (1/p, 1], 

The trace space 1/,^^ is then given by 

= {{9,d) G X : d,d = 0 on dn}, 

provided 

1 111 
p 2 p 2q 

otherwise one has to add 8,^9 = 0 in the definition of y^,^- In order to profit from the embedding 
(4.3) y^,^ -A cm X (IT; R"), 

we will always assume 

, ^ In 

1 > M > - + TT- 

p 2q 

Then by means of the assumptions stated before we obtain the following result. 

Theorem 2. Assume (R), Xjp nl2q < p < 1, and suppose that wg G Then the differential 

operator ^ 12 (^ 0 ) defined by A 2 (wo)w := A{wo,'V)w with domain D(yl 2 (wo)) := W has maximal 
Lp-regularity in Yg and thus also maximal Lp^p^-regularity in Yg. 


Proof. The proof is based on the results and techniques developed by Denk, Hieber, Priiss in [7] 
for the case /i = 1. By the results due to Priiss and Simonett EZ], these results extend to general 
fjL G (l/p + n/2q, 1], as the coefficients have enough regularity by the embedding (14.31) . 

a) The case fl = R" with constant coefficients 

In the sequel, we denote the covariable for thy z and that for x by The symbol A{£,) of A{wg, V) 
reads as 

aolCP + ai|c(0|2 -■j&o(Ao|^P+ 5^Ao|c(^)p)c(0'^ '' 

^blc{0 ^ieP + ^c(0®c(e) 


A0 = 


70 ' ^ ' 70 

where c(^) = ^ • Vdo- It is convenient to reduce this symbol for the variable Wred = [9,dred]^ where 
dred = c(^) ■ d. The reduced symbol ArediO becomes 

aolCP + ai|c(0P -*^o(Ao|CP + 9.Ao|c(e)n ^ 


Ared{f,) — 


70 70 

This symbol is homogeneous of second order and not strongly elliptic. However, it is normally elliptic 
in the sense of [7] as its spectrum satisfies ^{Aredif,)) C (0, 00 ) for each ^ 7 ^ 0. The latter can be seen 
by considering 

det(z + ArediO) = (z-\- Ool^p + ai|c(^)n(z + —+ ^^|c(0n 

7o 7o 


-bgbM0\H-\e 

70 


A0|^i2 , drXg 


7o 


c(e)n, 


which has two negative zeros, as 0 , 7 , k > 0 and A + 2 t 9,-A > 0. Therefore, by Section 6 of [7], the 
Lp-realization Ared of Ared has maximal Lp-regularity. This shows that whenever fg G Lp{J; Lg(R")) 
and fd G Lp(J; i7/(H; R")) are given, there is a unique solution 


of 


Wred = [9,dy Goi7p(J;i,(K”;R7)nLp(J;Lr,7R'*;R7) 

8tWred ~h Aj.^d'dJred — fred^ t 7 0, Wredih) = 0, 
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where fred = [fe, —*c(V) • fd]^- To obtain d, it remains to solve the problem 

dtd - ^ Ad = /] := fd + i^c{V)dred - 6ic(V)0, t > 0, d(0) = 0, 

70 70 

with maximal Lp-regularity of —A to obtain a unique solution 

d € oHliJ; K")) n Lp( J; K")), 

as A ,7 > 0 and G Lp(J; R")). This proves Theorem 2 in the case 17 = M" with constant 

coefficients. As detailed in Section 6 of |7], this assertion extends by perturbation and localization to 
variable coefficients, still in the case 17 = R". 


b) The case 17 = R" with constant coefficients 

It is convenient to replace x € R" by (x, y) € R"“^ x R+. On the symbolic level we have to replace 
^ by ^ — ivdy, where u denotes the outer normal at a boundary point of 17, and ^ • u = 0. Then c(^) 
becomes 

c(^ — ivdy) = (^ — iiydy) ■ Vdo = C ' • ^dody = ^ ■ Vdp, 

as 9iydo = 0. Therefore the symbol 7l(^) from Step a) is replaced by 

—Edy + z + 7l(^), 

where 

-iboc{Cy 
(Ao7o)/ 

Considering again the reduced variables Wred = (^, dred) with dred = c(^) • d, the reduced symbol 
becomes 


E = 


ao 

0 


~Ereddy + Z A.red{V)^ 


with 


E = 


ao 

0 


-ibo 

Ao7o 


To apply the half-space theory for normally elliptic operators in Section 7 of [7], we need to verify 
the corresponding Lopatinskii-Shapiro condition (LS) which states the following. 


If w G C'o(R+;R^) satisfies 

(LS) - Ereddyw{y) + {z + Ared{0)w{y) = 0, 2/> 0, apii;(0)=0, 

then w = 0. 


Condition (LS) can be proved without much pain. In fact, we observe that E^^y^{z -f ArediO) 
eigenvalues —< 0; otherwise z G C \ (— 00 ,0] would be a solution of 

0 = det(z -I- Ared{Cj + oj'^E red ) — det(z -|- Aredif + 

which by Step a) is impossible. Therefore B = {E~^^{z+Ared{£,))Y^'^ well-defined and has spectrum 
in C+. Thus, Wred{y) '■= e~^ywb is the unique stable solution of 

EredHy'^rediy) “t“ (.2^ T Ared{,^))'d^rediy) — 0; 2/ ^ b; '^redi^) — nJb- 

The Neumann condition implies 

0 = dyWrediff) — Bwb^ 

hence Wb = 0, as B is invertible, for all (z, ^ (0, 0), z G C \ (— oo, 0], ^ G M”, ( ■ iz = 0. Therefore, 

the techniques of Section 7 in |7] apply and show that the reduced problem has maximal Lp-regularity 
in R!(. 

As a result, given fe G Lp( J; Lq(R!J:)), fd G Lp( J; id,J(R!Ji; R”)), with fred = [fe, -*c(Va;) •/d]^, we 
find a unique solution Wred = [d, dred]^ of the problem 

Ht’dJred E Ared'^red — fred lu ; 'dJrediA) ~ 







16 


MATTHIAS HIEBER AND JAN TRUSS 


within the class 

Wred G oHliJ; L,(R;;R2)) nLp(J;ij2(K-;K2)), 

where Ared denotes the realization of Ared in with Neumann boundary condition. 

Next, we solve the remaining problem for d 

dtd - —Ad = fi in R", d(0) = 0, 

70 

d^d = 0 on 9K" 

with 

fd ■■= fd + i^ci\/)dred - bic{V)e, 

70 

in a similar way as in Step a) by employing maximal Lp-regularity for —A. This yields a unique 
solution 

d G oHliJ-, LgiRl; R”)) n Lp( J; ; R")). 

As S Lp( J; i7^(R” ; R")) we may differentiate the equation for d tangentially to obtain also 

V,dGoi^p(J; Lg(R!^;R(”-^)^”))nLp(J;i72(R-;R(—i)x-)). 

On the other hand, we may also take the derivative with respect to the normal variable y in order 
to obtain a problem with Dirichlet boundary conditions for v := dyd. We solve this with maximal 
Lp-regularity to obtain 

This proves Theorem 2 for the case 0 = R" with constant coefficients. As described in [7], Section 
7, this assertion extends by perturbation and localization to variable coefficients, still in the case 

n = R!;. 

c) General domains and variable coefhcients 

Here we follow the line given in [7], Section 8. We may use a perturbation argument to extend 
the result for the half-space to a bent half-space and then employ the localization method to prove 
Theorem 2 for general domains with -boundary. □ 


3. Local-Wellposedess 

We rewrite the above problem as an abstract quasilinear evolution equation of the form 
(4.4) z + A{z)z = F{z), t > 0, z{0) = zo- 

Here z = (u, w) = {u, 9, d) and we apply the Helmholtz projection P to the equation for u. The base 
space will be Xq := x Yq, where the subscript tr means solenoidal. Then with the generalized 

Stokes operator Ai(w) = —Pp(0,r)A, we define the regularity space by 

Ai := D(Ai) X Yi, D(Ai) = {u G R") n : m = 0 on 5H}. 

The operator A(z) is defined by A(z) = diag(Ai(i(;), A 2 (w)), and F{z) collects all lower order terms. 

In order to prove local well-posedness of (14.4L we may now resort to abstract theory, e.g. to the 
results by Kohne, Priiss, Wilke in and by LeCrone, Priiss and Wilke in |23| . 

Then, by Theorem 2 and by the maximal regularity of the generalized Stokes operator, see e.g. 
Bothe and Priiss P, A{z) has maximal Lp-regularity. For an interval J = [0,a], the solution space 
Ey(J) will be 

^ X,). 

The time-trace space X~^^y, of Ep( J) is given by 

= {UG n = 0} x 
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it satisfies 

^ X -A c(n)'*+i X c\n)^, 

provided 

(«) ; + 

Here denote as usual the Besov spaces; see e.g. Triebel [39]. Then A, F satisfy the requirements in 
the paper by LeCrone, Priiss and Wilke |23|, and so we have local well-posedness. If ^ + ^ + 5 < M < 1, 
the conditions of Kohne, Priiss, and Wilke |25| also hold. In particular, defining the state manifold 
of gH) by 

SM = {(u, e,d)€X^: 6» > 0, |d|2 = 1}, := X^^i, 

then SA4 is locally positive invariant for the semi-flow, total energy E is preserved, and the negative 
total entropy — N is a strict Lyapunov functional for the semi-flow on SA4 . Summarizing we have the 
following result 

Theorem 3. Assume (R), let p,q,p. be subject to (14.51) . and let zq G Xj^^. Then for some a = 
a(zo) > 0, there is a unique solution 

^ ^ ^ 0 ) n Xi), J = [0, a], 

of (14.4|) . i.e. (14.11) on J. Moreover, 

z GC{[0,a]■X^^p)f^C{{0,a]■X^), 

i.e. the solution regularizes instantly in time. R depends continuously on zq and exists on a maximal 
time interval J{zq) = [0,T'"(zo))- Moreover, 

t[^]zGLIpV(J;Xo)nLp,^(J;Xi), 

and \d(t,x )\2 = 1, E(t) = Eq, and —N is a strict Lyapunov functional. Furthermore, the problem 
(14.4p generates a local semi-flow in its natural state manifold SA4. 

4. The Generalized Principle of Linearized Stability 

Consider the autonomous quasilinear problem 

(4.6) z(t} + A(z(t))z(t) = F(z(t}}, t > 0, z(0) = Zq. 

Here we assume 

(4.7) (A,F) G C\V,B(Xi,Xo) x Xq), 

where V C Xj is open. Let £ C V n Xi denote the set of equilibrium solutions of (14.61) . which means 
that 

z G £ if and only if z G V n Xi, A(z)z = F(z). 

Given an element 2 * G £, we assume that z* is contained in an m-dimensional manifold of equilibria. 
This means that there is an open subset U C M™, 0 G U, and a C^-function f : U ^ Xi, such that 

(4.8) m)c£, m = ^*, Am)mo = Fmo), c&u. 

and the rank of 'tjj'{0) equals m. 

Let Aq denote the linearization of A{z)z — F{z) at z„, i.e. 

A[)h = A{zf}h F [^'(z*)h]z* — F'{zf)h. 

We call z* G f normally stable if the following conditions hold. 

(i) near z, the set f is a C^-manifold in Xi, dimf = m G No, 

(ii) the tangent space for £ at z, is isomorphic to N(Ho), 

(hi) 0 is a semi-simple eigenvalue of Aq, i.e. N(Ho) © R(7lo) = 

(iv) a(Ho) \ {0} C C+ = {C G C : ReC > 0}. 
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The following result is due to Priiss, Simonett and Zacher [38]. 

Theorem 4. Let \ <p < oo. Suppose z^, £ V Cl Xi is an equilibrium of and that [A, F) satisfy 
o and that A{z») has the property of maximal Lp-regularity. Assume further that z* is normally 
stable. 

Then z* is stable in Xj, and there exists <5 > 0 such that the unique solution z{t) of with 
initial value zq € X.y satisfying |zo — ^ exists on R+ and converges at an exponential rate in 

X.y to some Zoo & S as t ^ oo. 

It is worthwhile to note that in case m = 0, z* is necessarily isolated by (i). Then Theorem 4 reduces 
to the usual principal of linearized stability, as (ii), (iii), (iv) are equivalent to (j{Aq) C C+. 

5. Linear Stability of Equilibria 

The linearization of (|4.4I) . i.e. of (|4.1I) at an equilibrium z* = (0,0*, d*) is given by the operator 

A* = -diag((/i»/p)PA, (a*/pK*)A, (A»/ 7 *)A) 

in the base space Xq with domain D(A») = Xi. This operator has maximal Lp-regularity, it is 
the negative generator of a compact analytic Co-semigroup, and it has compact resolvent. So its 
spectrum consists only of countably many eigenvalues of finite multiplicity, which are all positive, 
hence stable, except for 0. The eigenvalue 0 is semi-simple, its eigenspace is given by 

N(A*) = {(0,d,d) : 0eR,deK"}, 

hence it coincides with the set of constant equilibria 6, when ignoring the constraint |d |2 = 1 and 
conservation of energy. Therefore each such equilibrium is normally stable. 

6. Nonlinear Stability 

We have stability with asymptotic phase for the equilibria of (14.41) . 

Theorem 5. Assume (R). Then any equilibrium z^: £ S of i4-4\ ) is stable in Xj. Moreover, for each 
Zt. £ £ there is e > Q such that if |zo — ij. — then the solution z of i4.4^ with initial value zq 

exists globally in time and converges at an exponential rate in X.y to some z^o G £■ 

This result is proved by means of the generalized principle of linearized stability, Theorem 4, above. 
In fact, by the previous section we know that each equilibrium z* = (0, 0*, d*) is normally stable. 

7. Long-Time Behaviour 

We conclude this paper with a result on the convergence of solutions to equilibria in the topology of 
the state manifold SA4. 

Theorem 6. Assume (R) and let z be the solution of i-e. of dH]), with initial value zq £ SA4. 

Then the following assertions hold. 

a) If we suppose 

sup [\z{t)\x^^,, + |l/0(t)|L„„] < oo, 
te(o,t+(zo)) 

then L'"(zo) = oo and z is a global solution. 

b) If z is a global solution, bounded in X^^p, and with 1/0 bounded, then z converges exponentially in 
SA4 to an equilibrium Zqo £ £ of |/./[ ), as t ^ oo. 

This result follows from abstract dynamical system arguments involving the strict Lyapunov func¬ 
tional — N, as well as the nonlinear stability result; see Kohne, Priiss, Wilke [22] ■ Note that, by a 
compactness argument, the converse of b) is also valid. 
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